DESCRIPTION 
FRACTAL STRUCTURE AND ITS FORMING METHOD 

This invention relates a fractal structure 
and its forming method, especially based on a novel 
principle. 

Background Art 

For application of a solid material to 
electronic or optical devices, physical properties of 
the material may restrict its applications. For 
example, in case of using a semiconductor material in a 
light -emitting device, it will be usable in a device of 
an emission wavelength corresponding to the band gap of 
the material, but some consideration will be necessary 
for changing the emission wavelength. Regarding 
physical properties related to semiconductor bands, 
controls by superlattices have been realized. More 
specifically, by changing the period of a superlattice , 
the bandwidth of its subband can be controlled to 
design an emission wavelength. 

Targeting on controlling many-electron-state 
structures by material designs, the Inventor proposed 
many-body effect engineering by quantum dot-based 
structures and has continued theoretical analyses ((1) 




U.S. Patent 5,430,309; (2) U.S. Patent 5,663,571; (3) 
U.S. Patent 5,719,407; (4) U.S. Patent 5.828,090; (5) 
U.S. Patent 5,831,294; (6) U.S. Patent 6.020,605; (7) 
J. Appl. Phys. 76, 2833(1994); (8) Phys . Rev. B51, 
10714(1995); (9) Phys. Rev. B51, 11136(1995); (10) J. 
Appl. Phys. 77, 5509(1995); (11) Phys. Rev. B53, 
6963(1996); (12) Phys. Rev. B53, 10141(1996); (13) 
Appl. Phys. Lett. 68, 2657(1996); (14) J. Appl. Phys. 
80, 3893(1996); (15) J. Phys. Soc . Jpn. 65, 3952(1996); 
(16) Jpn. J. Appl. Phys. 36, 638(1997); (17) J. Phys. 
Soc. Jpn. 66, 425(1997); (18) J. Appl. Phys. 81, 2693 
(1997); (19) Physica (Amsterdam) 229B, 146(1997); (20) 
Physica (Amsterdam) 237A, 220(1997); (21) Surf. Sci. 
375. 403(1997); (22) Physica (Amsterdam) 240B. 
116(1997); (23) Physica (Amsterdam) 240B. 128(1997); 
(24) Physica (Amsterdam) IE, 226(1997); (25) Phys. Rev. 
Lett. 80. 572(1998); (26) Jpn. J. Appl. Phys. 37, 
863(1998); (27) Physica (Amsterdam) 245B, 311(1998); 
(28) Physica (Amsterdam) 235B, 96(1998); (29) Phys. 
Rev. B59. 4952(1999); (30) Surf. Sci. 432. 1(1999); 
(31) International Journal of Modern Physics B. Vol. 
13, No. 21, 22, pp. 2689-2703, 1999). For example, 
realization of various correlated electronic systems is 
expected by adjusting a tunneling phenomenon between 
quantum dots and interaction between electrons in 
quantum dots. Let the tunneling transfer between 
adjacent quantum dots be written as t. Then, if 



quantum dots are aligned in form of a square lattice, 
the bandwidth of one electron state is T e ff = 4t . If 
quantum dots form a one-dimensional chain, the 
bandwidth of one electron state is T e ff=2t. In case of 
a three-dimensional quantum dot array, T e ff=6t. That 
is, if D is the dimension of a quantum dot array, the 
bandwidth of one electron state has been T e ff=2Dt. Here 
is made a review about half -filled (one electron per 
each quantum dot) Mott transition (also called Mott- 
Hubbard transition or Mott metal-insulator transition) . 
Let the effective interaction of electrons within a 
quantum dot be written as U e ff, then the Hubbard gap on 
the part of the Mott insulator is substantially 
described as A=U e ff-T e ff, and the Mott transition can be 
controlled by changing U e ff or t . As already proposed, 
the Mott-Hubbard transition can be controlled by 
adjusting U e ff or t, using a field effect, and it is 
applicable to field effect devices (Literatures (5), 
(6), (11) and (14) introduced above). 

On the other hand, reviewing the equation of 
A =Ueff-T ef f=Ueff-2Dt , it will be possible to control 
Mott-Hubbard transition by controlling the 
dimensionality D of the system. For this purpose, the 
Applicant already proposed a fractal-based structure 
that can continuously change the dimensionality, and 
have exhibited that Mott-Hubbard transition is 
controllable by changing the fractal dimensions. 



To enable designing of wider materials, it is 
desired to modify and control the dimension of 
materials by designing methods beyond the simple 
fractal nature. 

As physical systems in charge of information 
processing, intrinsic non-linearity is indispensable. 
As devices having been used for years, there are 
electronic devices using materials exhibiting non- 
linear responses to certain extents. For example, two- 
terminal devices exhibiting defferential negative 
resistances are one example of devices that are non- 
linear in current -voltage characteristics. Of course, 
MOS-FET, which is a three- terminal device, supports the 
modern technologies . By coupling these electronic 
devices having non- linear properties with a linear 
electronic circuit and constructing an information 
processing device having a non -linear property, any 
desired calculation can be executed. 

However, difficulties by high integration 
have become issues with such electronic circuits. 
Generation of heat is one of the problems, for example. 
The heat generation caused by intrinsic electric 
resistance is mandatory for producing non-linearity of 
an electronic device, and indispensable and essential 
for executing information processing. 

' To avoid this difficulty, trials have been 
made to decrease the number of devices by enhancing 



non-linearity of component devices. Progress of this 
scheme necessarily results in the need for component 
devices having non-linearity as strong as chaotic. 
When a chaotic classical system is quantized, what 
characterizes the behaviors of the quantum system is 
quantum chaos. 

On the other hand, as miniaturization of 
component devices progresses, electrons confined in a 
device result in behaving as quantum mechanical 
particles. Therefore, from this point of view, hopes 
are placed on components devices exhibiting quantum 
chaos. The Applicant has continued to theoretically 
demonstrate that, in a quantum system in a structure 
having a fractal configuration, the quantum chaos can 
be controlled by changing the fractal dimension 
characterizing the system. 

An object the invention intends to accomplish 
is to provide a fractal structure and its forming 
method capable of modulating and controlling the 
dimensionality of a material by a design method beyond 
the simple fractal nature. 

Another object the invention intends to 
accomplish is to provide a fractal structure and its 
forming method capable of controlling phase transition 
and chaos, in particular, quantum chaos, by a design 
method beyond the simple fractal nature. 
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The Inventor found, through concentrated 
researches toward solution of those issues, that a more 
complex fractal structure having a portion 
characterized by a plurality of fractal dimensions can 
be formed by changing the growth conditions with time 
during growth of the fractal structure. Especially in 
the process of growing a random fractal, it has been 
found that a distinctive nebula-like fractal-based 
complex structure made up of a mixture of a plurality 
of regions different in fractal dimension can be formed 
by changing growth conditions toward obtaining a higher 
fractal dimension after forming regions of a lower 
fractal dimension. Then, it has been found that, in 
fractal structures of this type, occurrence of phase 
transition, such as magnetic phase transition, and 
chaos, such as quantum chaos in an electron state, can 
be controlled. As a result of later detailed analysis, 
it has been found that there are fractal dimensions 
suitable for controlling these phenomena. 

The present invention has been made on the 
basis of those researches by the Inventor. 

That is, to overcome the above-indicated 
problems, according to the first aspect of the 
invention, there is provided a fractal structure having 
a plurality of regions different in fractal dimension 
characterizing the self -similarity , comprising: 



the fractal structure being grown from one or 
more origins under growth conditions providing a first 
fractal dimension in a portion of the growth process 
from the start point of time of growth to a first point 
and under growth conditions providing a second fractal 
dimension in another portion of the growth process from 
the first point of time to a second point of time. 

According to the second aspect of the 
invention, there is provided a fractal structure 
forming method for forming a fractal structure having a 
plurality of regions different in fractal dimension 
characterizing the self -similarity , comprising: 

growing the fractal structure from one or 
more origins under growth conditions providing a first 
fractal dimension in a portion of the growth process 
from the start point of time of growth to a first point 
and under growth conditions providing a second fractal 
dimension in another portion of the growth process from 
the first point of time to a second point of time. 

In the present invention , by adjusting the 
timing for changing the growth conditions from those 
for obtaining a first fractal dimension to those for 
obtaining a second fractal dimension, the nature of 
phase transition appearing in the fractal structure is 
controlled. For example, by adjusting the timing for 
changing the growth conditions, critical temperature 
for ferromagnetic phase transition appearing in the 



fractal structure is controlled. Further, by adjusting 
the timing for changing the growth conditions, the 
nature of chaos appearing in the fractal structure is 
controlled, or more particularly, quantum chaos in the 
electron state, for example, is controlled. The 
control of the quantum chaos in the electron state can 
be controlled with a high controllability by adjusting 
the timing for changing the growth conditions and 
additionally introducing a random magnetic field by 
addition of magnetic impurities. Typically, these 
regions forming the fractal structure exhibit a form of 
nebula as a whole. 

From the viewpoint of ensuring satisfactory 
control of the critical temperature of ferromagnetic 
transition appearing in the fractal structure and 
satisfactory control of the quantum chaos, or ensuring 
satisfactory control of the correlative electron 
system, the first fractal dimension D fi and the second 
fractal dimension D f2 are determined to be preferably 
Dfi>2.7 and D f2 <2.3, typically 2.7<D fl <;3 and lssD f2 <2.3, 
and more preferably 2 . 9^D f i^3 and l^D f2 <2.3. The upper 
limit value 3 of D fl corresponds to the dimension of a 
three-dimensional space whereas the lower limit value 1 
of D f2 is necessary for ensuring connectivity in the 
structure . 

In an embodiment that will be explained later 
in detail, growth conditions of the fractal structure 
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are represented by a of Equation (4). In actual 
growth, however, if the fractal structure is grown in 
the liquid phase, for example, natures of solvents 
employed for the growth are one of the growth 
conditions. That is, in this case, a plurality of 
regions different in fractal dimension can be formed by 
selecting appropriate solvents, respectively, in the 
growth process. 

According to the invention having the above - 
summarized configuration, a nebula-like fractal-based 
complex structure made up of a mixture of a plurality 
of regions different in fractal dimension from each 
other can be obtained by changing the growth condition 
of the fractal structure with time so as to first carry 
out the growth under growth conditions providing a low 
fractal dimension and to thereafter continue the growth 
under different growth conditions providing a higher 
fractal dimension. Then, in this fractal structure, 
the nature of phase transition occurring in the fractal 
structure can be controlled by, for example, adjustment 
of the timing for changing the growth conditions. 
Additionally, optimization of the fractal dimensions 
improves the controllability. 

Brief Description of Drawings 

Fig. 1 is a schematic diagram that shows a 
nebula-like fractal-based complex structure obtained by 



simulation according to the first embodiment of the 
invention; 

Fig. 2 is a schematic diagram that shows the 
nebula-like fractal-based complex structure obtained by 
simulation according to the first embodiment of the 
invention; 

Fig- 3 is a schematic diagram that shows the 
nebula- like fractal -based complex structure obtained by 
simulation according to the first embodiment of the 
invention; 

Fig, 4 is a schematic diagram that shows log- 
log plots of a distance r from the growth origin in the 
nebula-like fractal-based complex structure according 
to the first embodiment and the number of growth points 
N(r) contained in a ball with the radius r; 

Fig. 5 is a schematic diagram that shows 
spontaneous magnetization in the nebula-like fractal- 
based complex structure according to the first 
embodiment of the invention; 

Fig. 6 is a schematic diagram that shows a 
nebula- like fractal-based complex structure obtained by 
simulation according to the second embodiment of the 
invention; 

Fig. 7 is a schematic diagram that shows the 
nebula-like fractal-based complex structure obtained by 
simulation according to the second embodiment of the 
invention; 



Fig. 8 is a schematic diagram that shows the 
nebula-like fractal-based complex structure obtained by 
simulation according to the second embodiment of the 
invention; 

Fig. 9 is a schematic diagram that shows the 
nebula-like fractal-based complex structure obtained by 
simulation according to the second embodiment of the 
invention; 

Fig. 10 is a schematic diagram that shows 
quantum level statistics in the nebula-like fractal- 
based complex structure according to the second 
embodiment of the invention; 

Fig. 11 is a schematic diagram that shows 
quantum level statistics in the nebula- like fractal- 
based complex structure according to the second 
embodiment of the invention; 

Fig. 12 is a schematic diagram that shows a 
nebula- like fractal-based complex structure obtained by 
simulation according to the third embodiment of the 
invention; 

Fig. 13 is a schematic diagram that shows the 
nebula-like fractal-based complex structure obtained by 
simulation according to the third embodiment of the 
invention; 

Fig. 14 is a schematic diagram that shows the 
nebula-like fractal-based complex structure obtained by 
simulation according to the third embodiment of the 



invention; 

Fig. 15 is a schematic diagram that shows the 
nebula- like fractal-based complex structure obtained by 
simulation according to the third embodiment of the 
invention; 

Fig. 16 is a shematic diagram that shows 
quantum level statistics in the nebula-like fractal- 
based complex structure according to the third 
embodiment of the invention; 

Fig. 17 is a shematic diagram that shows 
quantum level statistics in the nebula-like fractal- 
based complex structure according to the third 
embodiment of the invention; 

Fig. 18 is a shematic diagram that shows a 
relation between a and the fractal dimension D f in a 
simplex fractal structure; 

Fig. 19 is a schematic diagram that shows 
spontaneous magnetization in a nebula-like fractal- 
based complex structure according to the fourth 
embodiment of the invention; 

Fig. 20 is a schematic diagram that shows 
spontaneous magnetization in the nebula-like fractal- 
based complex structure according to the fourth 
embodiment of the invention ; 

Fig. 21 is a schematic diagram that shows 
spontaneous magnetization in the nebula-like fractal- 
based complex structure according to the fourth 



embodiment of the invention; 

Fig. 2 2 is a schematic diagram that shows 
spontaneous magnetization in the nebula-like fractal- 
based complex structure according to the fourth 
embodiment of the invention; 

Fig. 2 3 is a schematic diagram that shows 
spontaneous magnetization in the nebula-like fractal- 
based complex structure according to the fourth 
embodiment of the invention; 

Fig. 24 is a schematic diagram that shows 
spontaneous magnetization in the nebula-like fractal- 
based complex structure according to the fourth 
embodiment of the invention; 

Fig. 2 5 is a schematic diagram that shows 
quantum level statistics in a nebula-like fractal-based 
complex structure according to the fifth embodiment of 
the invention; 

Fig. 26 is a schematic diagram that shows 
quantum level statistics in the nebula-like fractal- 
based complex structure according to the fifth 
embodiment of the invention; 

Fig. 2 7 is a schematic diagram that shows 
quantum level statistics in the nebula-like fractal- 
based complex structure according to the fifth 
embodiment of the invention; 

Fig. 28 is a schematic diagram that shows 
quantum level statistics in the nebula-like fractal- 



based complex structure according to the fifth 
embodiment of the invention ; 

Fig. 29 is a schematic diagram that shows 
quantum level statistics in the nebula- like fractal- 
based complex structure according to the fifth 
embodiment of the invention; 

Fig. 30 is a schematic diagram that shows 
quantum level statistics in the nebula- like fractal- 
based complex structure according to the fifth 
embodiment of the invention; 

Fig. 31 is a schematic diagram that shows 
quantum level statistics in the nebula-like fractal- 
based complex structure according to the fifth 
embodiment of the invention; 

Fig. 32 is a schematic diagram that shows 
quantum level statistics in the nebula- like fractal - 
based complex structure according to the fifth 
embodiment of the invention; 

Fig. 33 is a schematic diagram that shows a 
Berry-Robnik parameter p in the nebula- like fractal- 
based complex structure according to the fifth 
embodiment of the invention; 

Fig. 34 is a schematic diagram that shows a 
relation between energy and density of states in a 
nebula- like fractal-based complex structure according 
to the sixth embodiment of the invention; 

Fig. 35 is a schematic diagram that shows a 



relation between energy and density of states in the 
nebula-like fractal-based complex structure according 
to the sixth embodiment of the invention ; 

Fig. 36 is a schematic diagram that shows a 
relation between energy and density of states in the 
nebula-like fractal-based complex structure according 
to the sixth embodiment of the invention; 

Fig. 37 is a schematic diagram that shows a 
relation between energy and density of states in the 
nebula- like fractal-based complex structure according 
to the sixth embodiment of the invention; 

Fig. 38 is a schematic diagram that shows a 
relation between energy and density of states in the 
nebula- like fractal -based complex structure according 
to the sixth embodiment of the invention; 

Fig. 39 is a schematic diagram that shows a 
relation between energy and density of states in the 
nebula-like fractal-based complex structure according 
to the sixth embodiment of the invention; 

Fig. 40 is a schematic diagram that shows a 
relation between energy and density of states in the 
nebula-like fractal-based complex structure according 
to the sixth embodiment of the invention; 

Fig. 41 is a schematic diagram that shows a 
relation between energy and density of states in the 
nebula-like fractal-based complex structure according 
to the sixth embodiment of the invention; 



Fig. 42 is a schematic diagram that shows a 
relation between energy and density of states in the 
nebula-like fractal-based complex structure according 
to the sixth embodiment of the invention; 

Fig. 43 is a schematic diagram that shows a 
relation between energy and density of states in the 
nebula-like fractal-based complex structure according 
to the sixth embodiment of the invention; 

Fig. 44 is a schematic diagram that shows a 
relation between energy and density of states in the 
nebula-like fractal-based complex structure according 
to the sixth embodiment of the invention; 

Fig. 45 is a schematic diagram that shows a 
relation between energy and density of states in the 
nebula-like fractal-based complex structure according 
to the sixth embodiment of the invention; 

Fig. 46 is a schematic diagram that shows a 
relation between energy and density of states in the 
nebula- like fractal -based complex structure according 
to the sixth embodiment of the invention; 

Fig. 4 7 is a schematic diagram that shows a 
relation between energy and density of states in the 
nebula-like fractal-based complex structure according 
to the sixth embodiment of the invention; 

Fig. 48 is a schematic diagram that shows a 
relation between energy and density of states in the 
nebula-like fractal-based complex structure according 



to the sixth embodiment of the invention; 

Fig. 49 is a schematic diagram that shows a 
relation between energy and density of states in the 
nebula-like fractal-based complex structure according 
to the sixth embodiment of the invention; 

Fig. 50 is a schematic diagram that shows a 
relation between energy and density of states in the 
nebula-like fractal-based complex structure according 
to the sixth embodiment of the invention; and 

Fig. 51 is a schematic diagram that shows a 
relation between energy and density of states in the 
nebula-like fractal-based complex structure according 
to the sixth embodiment of the invention. 

BesJ^M ede for C aj,j. yiuy Otf t tT re ~ Tiiven Lloii 

Some embodiments of the invention will now be 
explained below. In the following embodiments, nebula- 
like fractal-based complex structures each made up of a 
mixture of portions characterized by a plurality of 
fractal dimensions are taken. These nebula-like 
fractal -based complex structures are formed by first 
forming a region having a lower fractal dimension and 
thereafter changing the growth conditions to provide a 
higher fractal dimension. 
First Embodiment 

(1) Formation of a nebula-like fractal-based complex 
structure 
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A method of forming a nebula-like fractal- 
based complex structure according to the first 
embodiment can be obtained by developing Dielectric 
Breakdown Model ((32) A. Erzan, L. Pietronero, A. 
Vespignani, Rev. Mod. Phys . 67, 545(1995); (33) L. 
Niemeyer, L. Pietronero, H.J. Wiesmann, Phys. Rev. 
Lett. 52, 1033(1984) (34) R. Ugajin, S. Hirata, and Y. 
Kuroki, Physica A, 278, 312(2000)). 

Here is defined a cubic lattice S in a three- 
dimensional space, and a scalar potential field (|>(ii, 
i 2 , i 3 ) is defined on its lattice site i 2 , ia)^S 

and called a potential. Let this potential obey the 
Laplace 1 s equation 

A0fe,/ 2 ,* 3 )-O (1) 
A pattern T n , which will be defined later, is a set of 
lattice sites on the three-dimensional lattice. T 0 
contains (0, 0, 0) alone, and T n+i is created by 
sequentially adding a single lattice site to T n 
according to the rule given below. 

Let the potential of each site contained in 
T n be 1 and let the potential at infinity be 0. That 
is, 

<f>(ii>h>h) = 0 when (h'h'h) 00 (2) 

^fc.'Wa)- 1 when (ii>h>h) GT n ( 3 > 
Under these boundary conditions, Equation (1) is solved 

to determine the potential. The lattice site, which 
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will be added to T n to construct T n+ i, is not contained 
in T n , and it is selected from the set of lattice 
sites, i.e. U n , nearest to T n . The number of lattice 
sites contained in U n is denoted by N n . 

Intensity of the electric fields for lattice 

sites (ii, m , i 2 ,nw ia.m) (where m=l, 2, N n ) in U n is 

defined as 

The probability that a certain site (ii, m , i 2 ,mr i3,m) in 
U n is selected is proportional to the electric field 
E ra ( a ) . That is, the probability is 

By repeating these operations, construction 
of T n is progressed. An ideal fractal will be a set of 
limits repeated infinitely as 

T„=limT n (6) 

n-*oo 

When a=l, the foregoing and the result of generation 
of the pattern by Diffusion limited aggregation ((35) 
T.A. Witten, Jr. and L.M. Sander, Phys . Rev. Lett. 47, 
1400(1984); Phys. Rev. B27, 5686(1983)). 

The nebula- like fractal-based complex 
structure according to the first embodiment is formed 
by changing the a parameter in accordance with the step 
n of the above-mentioned growth. That is, the above 
growth process is progressed using ai under l^n^x 1# a 2 



under Ti+l<sn5ST 2 , and a 3 under x 2 +l^n^T 3 . Simulation is 
performed below, taking fractal structure having two 
different fractal dimensions- Those characterized by 

cxi>a 2 are called nebula-like fractal-based complex 
structures, and simulation was carried out particularly 
for cases commonly fixing t 2 =10000, a x =2 and a 2 =0, and 
changing Xi- Results of the simulation are shown in 
Figs. 1, 2 and 3. Fig. 1 is for Ti=3000, Fig. 2 is for 
X!=5000, and Fig. 3 is for ^=7000. It will be 
understood from Figs. 1, 2 and 3 that regions of a low 
fractal dimension are formed while violently branching 
in an early period of growth, and thereafter, a layer 
of a higher fractal dimension grows on the branched 
structure. In the regions formed in the latter half of 
the growth, since the growth progresses entirely, it 
shows the tendency of smoothing the branched structure 
formed in the initial stage. This is just a nebula- 
like fractal-based complex structure. 

For the purpose of understanding the 
structure in greater detail, a calculation process of 
fractal dimensions is used. Let r represent the 
distance from the growth origin (0, 0, 0) and N(r) 
represent the number of growth points contained in a 
ball with the radius r. Then, using a as a factor of 
proportionality, if N(r) can be expressed as 

N(r)=ar Df (7) 
then D f is called the fractal dimension. Therefore, 



since logarithms of both sides are 

logN(r) = loga + D f logr ( 8 ) 

if log-log plots ride on a straight line, the structure 
can be regarded as a fractal, and its inclination is 
the fractal dimension. In Fig. 4, log-log plotting was 
made for various cases obtained by the above-mentioned 
growth experiment. It will be understood from Fig. 4 
that, under ^=3000, 5000 and 7000, the fractal 
dimension in the region with small log(r) is 
approximately D f =2.6, and the fractal dimension of the 
region with large log(r) is approximately D f =2.0. 
(2) Ferromagnetic phase transition in the nebula-like 
fractal-based complex structure 

A spin system describing a ferromagnetic 
structure is defined on the nebula-like fractal-based 
complex structure defined in (1). A review is made 
about the lattice point 

that is the origin of T n . In Equation (9), p=l, 2, 
n+1. By placing a spin on a lattice site r p ^T n , here 
is taken a spin system that can be described by the 
following Hamiltonian. 

(10) 

S p is the spin at the p site. A natural model of spin- 
to-spin interaction J p , q is: 



7 fl when |r,-r f |-l 
lO otherwise 



That is, spin-to-spin interaction exists only between 
nearest-neighbor sites. For the purpose of calculating 
spontaneous magnetization M at a finite temperature T, 
statistical mechanics of an equilibrium system is 
introduced, A partition function Z is defined as 
follows . 

Z = y e~ H/T (12) 
ft 

where {S p } in the symbol of the sum pertains to the sum 
regarding all spin states. Spontaneous magnetization 
is defined as the statistical average of spins as 
follows . 

-I n+l 

M=— — y(S„) (13) 
where the expected value <S P > is 

and n+l is the total number of spins. M is generally 
a vector quantity in a spin space, but its absolute 
value M = \M\ is calculated. 

Here is made a review about an Ising model. 
In an Ising model, only two states of 

5^=1^-1 (15) 
can exist. Let a mean field approximation be 
introduced into the Ising model. Spontaneous 



magnetization of the p-th site is written as \x p . In 
this system, since the molecule field varies with site, 
let it written as pL p . As an assumption of the mean 
field approximation, here is employed a molecule field 
that can be written by spontaneous magnetization of the 
nearest-neighbor site as 

Pp-^Jp*/** (16) 
This assumption simplifies the foregoing Hamiltonian to 

A self -consistent equation ensuring spontaneous 
magnetization obtained by using a partition function by 
the simplified Hamiltonian becomes results in 

^-tanh(^ p ) (18) 
and by numerically solving this equation, spontaneous 
magnetization of the system 

(19) 

y-o 

is obtained. 

Fig. 5 shows spontaneous magnetization in the 
nebula-like fractal-based complex structure shown in 
Figs. 1, 2 and 3. Note that the structure of Ti=0 is 
the same as a simple fractal structure having a single 
fractal dimension of a=0. Further, the structure of 
Ti=10000 is the same as a simple fractal structure 
having a single fractal dimension of a=2. It is seen 



that the critical temperature for spontaneous 
magnetization to disappear changes with x x . As apparent 
from Fig. 5, in this nebula-like fractal-based complex 
structure, the form of the magnetization curve does not 
change, and the critical temperature changes. It has 
been found, therefore, that various materials 
exhibiting various magnetic properties can be obtained 
by changing growth parameters in this nebula- like 
fractal-based complex structure, which is a fractal- 
based complex. 
Second Embodiment 

(1) Formation of a nebula-like fractal-based complex 
structure 

A method of forming the nebula-like fractal- 
based complex structure according to the second 
embodiment is the same as the method of forming the 
nebula-like fractal-based complex structure according 
to the first embodiment. Also in the second embodiment, 
assuming ai>a 2 , commonly fixing x 2 =10000, a x =2 and a 2 =0 # 
and changing Ti , simulation is carried out. Results of 
the simulation are shown in Figs. 6, 7, 8 and 9. Fig. 

6 is for Ti=2000, Fig. 7 is for Ti=4000, Fig. 8 is for 
Ti=6000, and Fig. 9 is for Xi=8000. As seen from Figs. 
6, 7, 8 and 9, similarly to the first embodiment, since 
regions of a low fractal dimension are formed while 
violently branching in an early period of growth; a 
layer of a higher fractal dimension thereafter grows on 



the branched structure; and the growth progresses 
entirely in the regions formed in the latter half of 
the growth, it shows the tendency of smoothing the 
branched structure formed in the initial stage. 
(2) Electron system on the nebula-like fractal-based 
complex structure 

A quantum system of a single particle is 
defined on the nebula- like fractal-based complex 
structure defined in (1). Assume a lattice site shown 
below, which is the origin of T n . 

r p -(i^i^J^GT. (20) 
where p=0, 1, 2, n. Here is defined an operator c* 

that creates a quantum at a lattice site r p ^T n . Of 
course, an anticommutative relation 



(21) 

is established. Here the Hamiltonian H is defined as 

" = -2'-% < 22 > 



p,q 



Here is employed as the transfer t 



1 when Yp ~ r q I = 1 



(23) 



0 otherwise 
In this model, hopping is possible only between 
nearest -neighbor sites . 

When £ m denotes the eigenenergy of the 
Hamiltonian H and I m) denotes the eigenvector, 

H\m) = e m \m) (24) 



where m=0 , 1, 2, n. 

First, n+l quantum levels £ m are 
standardized such that spacing between nearest -neighbor 
levels becomes 1 in average* That is, 

co j -e r e M (25) 
However, when j = l, 2, n, by using 

a7 = -Ya) y (26) 
"Pi 

it is converted into new levels 

e Q =0 (27) 



; {A J 

Here, 



^co^YQ; (28) 



CO 

Q,--zr (29) 

CO 



The density of states of the system is defined by 

(30) 

n +1 ^ 



m=0 

and the staircase function 



A(6)-f drypfr) (3D 



is calculated. The staircase function obtained is 
converted by a technique called "unfolding" such that 
the density of states becomes constant in average- By 
using quantum levels obtained in this manner, nearest - 
neighbor level spacing distribution P(s) and A 3 
statistics of Dyson and Metha are calculated as quantum 



level statistics. As taught in literatures ((36) 
(37) ), by using these statistics, it can be detected 
whether quantum chaos has been generated or not. It is 
also known that a quantum chaotic system is sensitive 
to perturbation from outside similarly to the classical 
chaotic system, and analysis of quantum chaos is 
important as a polestar of designs of non-linear 
materials . 

In case of an integrable system, nearest- 
neighbor level spacing distribution P(s) and A 3 
statistics are those of Poisson distribution 

P P (s) = e- s (32) 

A 3 («) = ^ (33) 

In case of a quantum chaotic system, it becomes GOE 
(Gaussian orthogonal ensemb) distribution 

Poo E (s) = fe-^ (34) 
A 3 (n) = ^[log(2™) + y-^-f| + o(«- 1 ) (35) 

Jt o 4 

where 7 is the Euler's constant. 

Since the nebula-like fractal-based complex 
structure analyzed here is one obtained by growth 
experiment of n=10000, this quantum system includes 
n+l=10001 eigenstates. Based on energy eigenvalues 
concerning 1501 states from the 510-th to 2001-th 
eigenstates from the ground state, among those 



eigenstates, the following quantum level statistics was 
calculated. Figs. 10 and 11 show quantum level 
statistics in a nebula-like fractal-based complex 
structure of {a lf a 2 ) = (0, 2). Fig. 10 shows P(s) and 
Fig. 11 shows A 3 statistics. In case of Ti=10000, the 
structure is the same as the simple fractal of a=2, and 
the fractal dimension is near 2(D f -2.16). Therefore, 
the system behaves as an integrable system. In this 
case, the quantum level statistics is that of Poisson 
distribution. On the other hand, in case of x x =0, the 
structure is the same as the simple fractal of a=0, and 
the fractal dimension is near 3(D f —2.91). Therefore, 
the system behaves as a quantum chaotic system, and the 
quantum level statistics is that of GOE distribution. 
As shown in Figs. 10 and 11, as Xi increases from 0 to 
10000, the quantum level statistics changes from GOE 
distribution to Poisson distribution. Therefore, wide 
kinds of quantum systems can be realized by setting Xi 
at predetermined values . 

Through a review of the changes in quantum 
level statistics in the nebula-like fractal-based 
complex structure, it is recognized that the quantum 
chaotic properties begin to decrease with a relatively 
small value of Xi. This suggests that the branching 
property still remains sufficiently in the complex 
structure obtained and the chaotic property of the 
quantum chaos is weakened. Taking it into 
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consideration, it is recognized that a plurality of 
regions conjoin in the nebula-like fractal-based 
complex structure and form a complex structure as a 
whole . 
5 Third Embodiment 

(1) Formation of a nebula-like fractal-based complex 
structure 

A method of forming the nebula-like fractal- 
based complex structure according to the third 
,^10 embodiment is the same as the method of forming the 

y nebula-like fractal-based complex structure according 

Q 

to the first embodiment. In the third embodiment, 
( -S assuming ai>a 2 , commonly fixing t 2 =10000, ai=2 and a 2 =0, 

m 

'"'"4 and changing Ti, simulation is carried out. Results of 

is 

Q15 the simulation are shown in Figs. 12, 13, 14 and 15. 

□ Fig. 12 is for t 1 =2000. Fig. 13 is for Ti=4000, Fig. 14 

ru 

rg is for Ti=6000, and Fig. 15 is for Ti=8000. As seen 

!y from Figs. 12, 13, 14 and 15, similarly to the first 

embodiment, since regions of a low fractal dimension 
20 are formed while violently branching in an early period 

of growth; a layer of a higher fractal dimension 
thereafter grows on the branched structure; and the 
growth progresses entirely in the regions formed in the 
latter half of the growth, it shows the tendency of 
2 5 smoothing the branched structure formed in the initial 

stage . 

(2) Electron system on the nebula-like fractal-based 
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complex structure 

A quantum system of a single particle is 
defined on the nebula- like fractal-based complex 
structure defined in (1). Assume a lattice site shown 
below, which is the origin of T n . 

r p =(i hp J 2 , p J 3tP )eT n (36) 
where p=l, 2, ••, n+1. Here is defined an operator c* 
that creates a quantum at a lattice site r p ^T n . Of 
course, an anticommutative relation 



(37) 

is established. Here the Hamiltonian H is defined as 

H =-yt BO c\c a (38) 



PA 



Here is employed as the transfer t p , q , 

1 when |r - rj = 1 



p i\ 
0 otherwise 



(39) 



where 9p f q=-e q , p is a random real number satisfying 



0<e pq <7jt (40) 
In this model, hopping is possible only between 
nearest -neighbor sites. Then, along with the hopping, 
phase factors 9 p , q , which are random from one site to 
another, are added. If the phase factor is integrated 
in the loop making one turn around a lattice point, a 
magnetic flux passing through the loop is obtained. 
This means a magnetic field is locally introduced to 
the random distribution of 0<9 p , q <2jt. This magnetic 



field is absolutely random in both intensity and 
direction, and in spatial average, it becomes a zero 
magnetic field and never breaks the fractal property of 
the system. 

5 When £ m denotes the eigenenergy of the 

Hamiltonian H and I m) denotes the eigenvector, 

H\m) = e m \m) (41) 

where m=0, 1, 2, n. 

First, n+1 quantum levels £ m are 
10 standardized such that spacing between nearest -neighbor 

levels becomes 1 in average. That is, 

co j -e r e hl (42) 
However, when j = l, 2, n, by using 

O d7 = ~2 (0 J (43) 

□ 

nil 5 it is converted into new levels 



□ 
CO 



ru 



<» ft ft 

Here , 



e 0 =0 (44) 

1 m m 

^-7=2^ -1 Q > < 45) 



CO, 

Q,=-d- (46) 



1 co 



20 The density of states of the system is defined by 



and the staircase function 



31 



A(e) = f £ drj Pin) (48) 

J — 00 

is calculated. The staircase function obtained is 
converted by a technique called "unfolding" such that 
the density of states becomes constant in average. By 
using quantum levels obtained in this manner, nearest - 
neighbor level spacing distribution P(s) and A 3 
statistics of Dyson and Metha are calculated as quantum 
level statistics. As taught in a literature ((36) and 
(37) shown above), by using these statistics, it can be 
detected whether quantum chaos has been generated or 
not. It is also known that a quantum chaotic system is 
sensitive to perturbation from outside similarly to the 
classical chaotic system, and analysis of quantum chaos 
is important as a polestar of designs of non-linear 
materials . 

In case of an integrable system, nearest- 
neighbor level spacing distribution P(s) and A 3 
statistics are those of Poisson distribution 

P P (s) = e~ s (49) 

A 3 (n) = ^ (50) 

In case of a quantum chaotic system under a magnetic 
field, it becomes GUE (Gaussian unitary ensemble) 
distribution 

Ws) = ^-e- 4W ' T (51) - 

Jt 



Lfl 



•sss 



log(2;m) + y - ^ 



+ 0(0 (52) 



where 7 is the Euler's constant. 

Since the nebula-like fractal-based complex 
structure analyzed here is one obtained by growth 
5 experiment of n=10000, this quantum system includes 

n+l=10001 eigenstates. Based on energy eigenvalues 
concerning 1501 states from the 510-th to 2001-th 
eigenstates from the ground state, among those 

M 

□ eigenstates, the following quantum level statistics was 

Q 

mlO calculated. Figs, 16 and 17 show quantum level 



statistics in a nebula-like fractal-based complex 
structure under a random magnetic field of (a i# a 2 ) = (2, 
0). Fig. 16 shows P(s) and Fig. 17 shows A 3 
statistics. In case of Ti=10000, the structure is the 



□ 

I'M 

f .~15 same as the simple fractal of a=2, and the fractal 

u 

TU dimension is near 2(D f -2.16). Therefore, the system 

behaves as an integrable system. In this case, the 
quantum level statistics is that of Poisson 
distribution. On the other hand, in case of Ti=0, the 

20 structure is the same as the simple fractal of a=0 , and 

the fractal dimension is near 3(D f -2.91). Therefore, 
the system behaves as a quantum chaotic system, and the 
quantum level statistics is that of GUE distribution. 
As shown in Figs. 16 and 17, as x 1 increases from 0 to 

25 10000, the quantum level statistics changes from GUE 

distribution to Poisson distribution . Therefore , wide 
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kinds of quantum systems can be realized by setting x x 
at predetermined values . 
Fourth Embodiment 

(1) Formation of a nebula-like fractal-based complex 
structure 

A method of forming the nebula-like fractal- 
based complex structure according to the fourth 
embodiment is the same as the method of forming the 
nebula-like fractal-based complex structure according 
to the first embodiment. As growth conditions, 
however, while fixing Ti=5000 and t 2 =10000, (cti, oc 2 ) was 
changed variously within the range satisfying ai>a 2 for 
carrying out numerical experiment of growth. 

Regarding methods of forming a single fractal 
structure by a dielectric breakdown model, it is known 
that changes of a invite changes of the fractal 
dimension D f of the fractal structure to be formed. 
Fractal dimensions having obtained by simulation are 
shown in Fig. 18 (see Literature (34)). As apparent 
from Fig. 18, as a increases, D f decreases. When 
a<0.5, D f >2.7 is obtained. When a>l, D f <2.3 is 
obtained. 

(2) Ferromagnetic phase transition in a nebula-like 
fractal-based complex structure 

A spin system describing a ferromagnetic 
structure is defined on the nebula-like fractal-based 
complex structure defined in ( 1 ) . A review is made 



about the lattice point 

r p - 0'^.«2 iP . f 3.p) e7 i. (53) 
that is the origin of T n . In Equation 53, p=0, 1, 2, 

n. By placing a spin on a lattice site r p ^T n , here 
is taken a spin system that can be described by the 
following Hamiltonian . 



H = ~yj Da S -S a (54) 

/ j PA P <? x ' 



PA 



S p is the spin at the p site. A natural model of spin- 
spin interaction J p , q is: 



" 3 10 J = J 



1 when r\, - r J = 1 , 

(55) 



P q 

0 otherwise 



|»0 That is, spin-spin interaction exists only between 

SI nearest -neighbor sites. For the purpose of calculating 

q spontaneous magnetization M at a finite temperature T, 

S3 

i:3 statistical mechanics of an equilibrium system is 

ro 

r -^15 introduced. A partition function Z is defined as 

fy follows. 

Z = V e- H/T (56) 
(?) 

where {S p } in the symbol of the sum pertains to the sum 
regarding all spin states. Spontaneous magnetization 
20 is defined as the statistical average of spins as 

follows. 



where the expected value <S P > is 
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( s >)-v2 s > e ~" T (58) 



and n+1 is the total number of spins. M is generally 
a vector quantity in a spin space, but its absolute 
value Af =|A/| is calculated. 

Here is made a review about an I sing model. 
In an Ising model, only two states of 

S p =lor -1 (59) 
can exist. Let a mean field approximation be 
introduced into the Ising model. Spontaneous 
magnetization of the p-th site is written as (jt p . In 
this system, since the molecule field varies with site, 
let it written as Ji . As an assumption of the mean 
field approximation, here is employed a molecule field 
that can be written by spontaneous magnetization of the 
nearest-neighbor site as 

^-S'm". (60) 

This assumption simplifies the foregoing Hamiltonian to 

H MF =-%Jl p o p (61) 

A self -consistent equation ensuring spontaneous 
magnetization obtained by using a partition function by 
the simplified Hamiltonian becomes \x p results in 

^ - tanh(0F,) (62) 
and by numerically solving this equation, spontaneous 
magnetization of the system 
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M lsin g =^-r5X (63) 

n + 1 ^ 

y=o 

is obtained. 

Fig, 19 shows spontaneous magnetization in 
case of a 2 =0 fixed and cti changed. Here, ax=0 
corresponds to spontaneous magnetization in the same 

structure as the simple fractal of a=0. As shown in 
Fig. 19, as cti increases, the phase transition 
temperature of the nebula-like fractal-based complex 
structure is remarkably changed. 

Fig. 20 shows spontaneous magnetization in 
case of cx 2 =0.2 fixed and a x changed. Here, cti=0.2 
corresponds to spontaneous magnetization in the same 
structure as the simple fractal of a=0.2. The 
modulatory property is smaller as compared with Fig. 
19, but still remains significantly large. 

Fig. 21 shows spontaneous magnetization in 
case of ct 2 =0.4 fixed and a x changed. Here, ai=0.4 
corresponds to spontaneous magnetization in the same 
structure as the simple fractal of a=0.4. As seen from 
Fig. 21, the modulatory property is remarkably smaller 
as compared with Figs. 19 and 20. 

Fig. 22 shows spontaneous magnetization in 
case of ai=2 fixed and <x 2 changed. Here, a 2 =2 
corresponds to spontaneous magnetization in the same 
structure as the simple fractal of a=2. Fig. 23 shows 
spontaneous magnetization in case of ai = l fixed and a 2 



changed. Here, a 2 =l corresponds to spontaneous 
magnetization in the same structure as the simple 

fractal of a=l. Fig. 24 shows spontaneous 
magnetization in case of ai=0.6 fixed and a 2 changed. 
Here, a 2 =0.6 corresponds to spontaneous magnetization 
in the same structure as the simple fractal of a=0.6. 
As seen from Figs. 22, 23 and 24, the phase transition 
temperatures are controlled by a 2 to a certain extent. 

Generalizing the above aspects of spontaneous 
magnetization, when a nebula-like fractal-based complex 
structure satisfies the condition of approximately 
D f <2.3 in the growth process of lssnssri, i.e. in the 
initial stage of the growth, and the growth condition 
of approximately D f >2.7 in the growth process of 
Ti+l<sn£T 2 , i.e. in the final stage of the growth, the 
ferromagnetic phase transition temperature in the 
nebula-like fractal-based complex structure will be 
controlled well. 
Fifth Embodiment 

(1) Formation of a nebula-like fractal-based complex 
structure 

A method of forming the nebula- like fractal- 
based complex structure according to the fifth 
embodiment is the same as the method of forming the 
nebula-like fractal-based complex structure according 
to the first embodiment. As growth conditions, 
however, fixing Ti=5000 and x 2 =10000 and using various 
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combinations of (<Xi, a 2 ) satisfying ai>a 2 , particularly- 
using 0, 0.2, 0.4, 0.6, 1 and 2 as cti and a 2 , numerical 
experiment of growth was carried out. 

(2) Electron system on the nebula-like fractal-based 
complex structure 

In the same process as demonstrated by 
Equations (20) through (35) of the second embodiment, a 
quantum system of one particle is defined on the 
nebula-like fractal-based complex structure defined in 
(1) • 

Since the nebula-like fractal-based complex 
structure analyzed here is one obtained by growth 
experiment of n=10000 # this quantum system includes 
n+l=10001 eigenstates. Based on energy eigenvalues 
concerning 1501 states from the 510-th to 2001-th 
eigenstates from the ground state, among those 
eigenstates, the following quantum level statistics was 
calculated. Figs. 2 5 and 26 show quantum level 
statistics in a nebula-like fractal-based complex 
structure under (cxi, a 2 ) = (x, 2) where x=0, 0.2, 0.4, 
0.6, 1 and 2. Fig. 25 shows P(s) and Fig. 26 shows A 3 
statistics. In case of (Xi=0, the structure is the same 
as the simple fractal of a=0, and the fractal dimension 
is near 3(D f -2.91). Therefore, the system behaves as 
a quantum chaotic system. In this case, the quantum 
level statistics is that of GOE distribution. As a x 
increases, the quantum level statistics changes away 
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from that of GOE distribution and approaches that of 
Poisson distribution. However, even after reading 
cti=2, a large difference from Poisson distribution 
still remains - 

Fig. 27 shows A 3 statistics in a nebula-like 
fractal-based complex structure under (cx lf a 2 )=(x, 0.2) 
where x=0.2, 0.4, 0.6, 1 and 2. In case of cti=0.2, the 
structure is the same as the simple fractal of cx=0.2, 
and behaves as a quantum chaotic system. As ai 
increases, the quantum level statistics changes away 
from that of GOE distribution and approaches that of 
Poisson distribution. Figs. 28 shows A 3 statistics in 
a nebula- like fractal -based complex structure under 
(cti# a 2 ) = (x, 0.4) where x=0.4, 0.6, 1 and 2. In case 
of ai=0.4, the structure is the same as the simple 
fractal of a=0.4, and it is already apart from that of 
GOE distribution too much to be a quantum chaotic 
system. Here again, as cti increases, the quantum level 
statistics changes away from that of GOE distribution 
and approaches that of Poisson distribution. 

Figs. 29 and 30 show quantum level statistics 
in a nebula-like fractal-based complex structure under 
(oti, a 2 ) = (2, x) where x=0, 0.2, 0.4, 0.6, 1 and 2. 
Fig. 29 shows P(s) and Fig. 30 shows A 3 statistics. In 
case of a 2 =0, it is the same as ai=2 of Figs. 25 and 
26. In case of a 2 =2, the structure is the same as the 
simple fractal of a=2, and the fractal dimension is 



near 2(D f -2.16). Therefore, the system behaves as an 
integrable system. As a 2 decreases, the quantum level 
statistics changes away from that of Poisson 
distribution and approaches GOE distribution. 

Fig. 31 shows A 3 statistics in a nebula-like 
fractal-based complex structure under (a i# a 2 )=(l, x) 
where x=0 , 0.2, 0.4, 0.6 and 1 . Fig . 32 shows A 3 
statistics in a nebula-like fractal-based complex 
structure under (a i# a 2 )=(0.6, x) where x=0, 0.2, 0.4 
and 0.6. As a 2 decreases, the quantum level statistics 
changes away from that of Poisson distribution and 
approaches that of GOE distribution. 

For quantitative evaluation of the above - 
reviewed controllability, the Berry-Robnik parameter p 
is used ((38) M.V. Berry and M. Robnik, J. Phys . A 
(Math. Gen.) 17, 2413(1984)). First, when p=l-p, 

P 2 (s,p) = P 2 e-erf(^) - { 2 pp)^y-^ ( 64 ) 

is introduced , where 

er f{x) = ~r x dre- Tl (65) 
-J it 

was used. This function P 2 (s, p) coincides with P(s) 
of Poisson distribution under p=l, and coincides with 
P(s) of GOE distribution under p=0. That is, by 
changing p from 0 to 1 , quantum level statistics from 
those of quantum chaotic systems to those of integrable 
systems can be interpolated. The Berry-Robnik 



parameter is the value of p obtained by optimum 
approximation of P(s) obtained y numerical calculation 
with P 2 (s, p) shown above. Within the range of 
semiclassical approximation, p is the ratio in volume 
of regular regions (integrable systems and regions 
capable of perturbing development therefrom) in a phase 
space - 

Fig. 33 shows the Berry-Robnik parameter p in 
nebula- like fractal-based complex structures. (cti, 
a 2 )=(x / 0) is the Berry-Robnik parameter in the case 
where cti is placed on the abscissa while fixing a 2 =2. 
It is apparent from Fig. 33 that widely various quantum 
systems from quantum chaotic systems to integrable 
systems can be realized by setting {a lf a 2 ) to 
predetermined values. 
Sixth Embodiment 

(1) Formation of a nebula-like fractal-based complex 
structure 

A method of forming the nebula- like fractal - 
based complex structure according to the sixth 
embodiment is the same as the method of forming the 
nebula- like fractal-based complex structure according 
to the first embodiment. As growth conditions, 
however, fixing Ti=5000 and x 2 =10000 and using various 
combinations of (cti, ot 2 ) satisfying a x >a 2# particularly 
using 0, 0.2, 0.4, 0.6, 1 and 2 as a x and a 2 , numerical 
experiment of growth was carried out . 



+ 



(2) Correlative electron system on the nebula-like 
fractal-based complex structure 

An electron system is defined on the nebula- 
like fractal-based complex structure defined in (1). A 
review is made about the lattice point 

that is the origin of T n . In Equation (66), p=0, 1, 2, 
n. An operator c* CT for generating an electron of a 
spin a is defined at the lattice point r p E.T n . Of 
course, there is the anticommutative relation 

Here is defined a single-band Hubbard Hamiltonian H of 
the electron system as follows . 



H=tVXJlc. a+U yn j ,n. i (68) 



i,J>o J 

Assuming that electrons are movable only 
among nearest-neighbor sites, the following is employed 

aS ^p,q- 

1 when |r — r I = 1 



P>9 



P <7 I 

0 otherwise 



(69) 



Additionally, the spin o electron density operator of 
the j-th site, n ja =c) a c ja , and their sum, n } =^ a ^j, a • 
are defined. 

For the purpose of defining a temperature 
Green's function, here is introduced a grand canonical 
Hamiltonian K = H - uN where N = ^V« y ' In the half ~ 
filled taken here, chemical potential is M=U/2. The 
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half -filled grand canonical Hamiltonian can be 
expressed as 



Operators t jJa , j jJ[o , u i and d io are defined beforehand 
as 



^. c. + c T 



'/,.a =CJ ><7 C,. 0 +C7 (£J ^, CT (71) 



7Mo = C ,,a C i,a -c'aC/.c (72) 

«. = KAA.A\ + £ i.A.AA.x {73) 

If the temperature Green function is defined for 

operators ^4 and 5 given, taking r as imaginary time, 
it is as follows. 

(a ; b) = -f fi Q dr(T r A(r)By^ (75) 

The on-site Green function 

G ha {i<o n )-(c ha ;cl a ) (76) 
is especially important. 

Imaginary time development of the system is 
obtained by the Heisenberg equation 

-4-A(t) = [k,a] (77) 
ax 

As the equation of motion of the on-site Green 
function , 



is obtained. Then, the approximation shown below is 



introduced, following Gros ((35) C. Gros, Phys . Rev. 
B50, 7295(1994)). If the site p is the nearest- 
neighbor site of the site j , the resolution 

^t(c p y,c^)(c } y,cl) ( 79 ) 

is introduced as the approximation. This is said to be 
exact in case of infinite-dimensional Bethe lattices, 
but in this case, it is only within approximation. 
Under the approximation, the following equation is 
obtained. 

(io>n -tX,o)Gj, a =l + y(rf y ,. a c ; , a ;c; a ) ( 80 ) 

where 

r y,*=2 A - G - (81) 

P 

was introduced. To solve the equation, obtained, 

{^j,-a^j,o'^) t a) has to be analyzed. In case of half-filled 

models, this equation of motion is 

ico Id. c :c f ) = —G- (1 • c. :c t \ 

p 

+ i l\j(d,-oC p , 0 -Ao) (82) 
Here again, with reference to the Gros theory, 

approximation is introduced. It is the following 

translation. 

(hj.-o^l) - -<G p ,_ 0 (d._J J , 0 ;cl) ( 83 ) 



(dj.-Jj* ; c[ a ) - tG„ (d L _ a c ji0 ; c J, a ) (84) 



By executing this translation, the following closed 
equation is obtained. 

(iC0 n - 1%„ )G jt0 = 1 + — -|^L_ G . o ( 85 ) 

Here is assumed that there is no dependency on spin, 
that is, assuming 

the following calculation is executed. This is 
because, when analytic continuation is conducted as 
ico n -^co+i 5 for small 6, 

Dj(co) = -lmGj(a) + id) (87) 
becomes the local density of states of the site j , and 

D(oj) = _J_y D (co) (88) 

n +1 *f J 

becomes the density of states of the system. For later 
numerical calculation of densities of states, 6=0.0001 
will be used. 

Regarding the nebula-like fractal-based 
complex structure obtained in (1), fixing t=l, Ti=5000 
and x 2 =10000, densities of states are calculated. Figs. 
34, 35 and 36 show densities of states of half -filled 
electron systems (10001 electrons) in the nebula-like 
fractal-based complex structure under (a lt a 2 )=(x, 0) 
where x=0, 0.2, 0.4, 0.6, 1 and 2. U=6 is used in Fig. 
34, U=7 in Fig. 35, and U=8 in Fig. 36. In case of 



(*!=(), the structure is the same as the simple fractal 
of a=0, and the fractal dimension is near 3(D f -2.91). 
In Fig. 34 , for all values of a 1# densities of states 
D(0) under co=0 are finite, and the electron system is 
in a metal phase. As cti increases, D(0) gradually 
decreases, and the electron system approaches an 

insulator. In Fig. 35, when a x <l, the state of density 
D(0) under co=0 is finite, and the electron system is 
in a metal phase. However, when a x =2 , D(0) disappears 
substantially, and Mott insulator transformation 
occurs. In Fig. 36, as cii increases, disappearance of 
D(0) occurs, that is, typical Mott transition occurs. 

Figs. 37, 38 and 39 show densities of states 
of half -filled electron systems (10001 electrons) in 
the nebula- like fractal-based complex structure under 
a 2 )=(x, 0.2) where x=0.2, 0.4, 0.6, 1 and 2. U=6 
is used in Fig. 37, U=7 in Fig. 38, and U=8 in Fig. 39. 
In case of a x =0 . 2 , the structure is the same as the 
simple fractal of ct=0.2. In Fig. 37, for all values of 
cti, densities of states D(0) under co=0 are finite, and 
the electron system is in a metal phase. As cx x 
increases, D(0) gradually decreases, and the electron 
system approaches an insulator. In Fig. 38, when <Xi<l, 
the state of density D(0) under co=0 is finite, and the 
electron system is in a metal phase. However, when 
cti=2, D(0) disappears substantially, and Mott insulator 
transformation occurs. In Fig. 39, it is recognized 



the system behaves as a Mott insulator under all values 
of ai, but its insulation performance, i.e. electron- 
to-electron correlative effect, is enhanced as cti 
increases. This is because width of the effective 
Hubbard gap increases along with the increase of cti . 

Figs. 40, 41 and 42 show densities of states 
of half -filled electron systems (10001 electrons) in 
the nebula- like fractal-based complex structure under 
(cxi, a 2 ) = (x, 0.4) where x=0.4, 0.6, 1 and 2. U=6 is 
used in Fig. 40, U=7 in Fig. 41, and U=8 in Fig. 42. 
In case of ai=0.4, the structure is the same as the 
simple fractal of a=0.4. In Fig. 40, for all values of 
cti, densities of states D(0) under co=0 are finite, and 
the electron system is in a metal phase. As a x 
increases, D(0) gradually decreases, and the electron 
system approaches an insulator. In Fig. 41, as cti 
increases, disappearance of D(0) occurs, that is, 
typical Mott transition occurs. In Fig. 42, for all 
values of a lt the systems behave as Mott insulators, 
but it is recognized that the insulation performance, 
i.e. electron-to-electron correlative effect, increases 
along with the increase of (Xi. This is because the 
width of the effective Hubbard gap becomes larger as ai 
increases . 

Figs. 43, 44 and 45 show densities of states 
of half -filled electron systems (10001 electrons) in 
the nebula-like fractal-based complex structure under 
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(ct 1# a 2 ) = (2, x) where x=0, 0.2, 0.4, 0.6, 1 and 2. U=6 
is used in Fig. 43, U=7 in Fig. 44, and U=8 in Fig. 45. 
In case of a 2 =2, the structure is the same as the 
simple fractal of a=2, and the fractal dimension is 
near 2(D f —2.16). In Fig. 43, for all values of a 2 , 
densities of states D(0) under co=0 are finite, and the 
electron system is in a metal phase. As a 2 increases, 
D(0) gradually decreases, and the electron system 
approaches an insulator. In Figs. 44 and 45, for all 
values of a 2 , the systems behave as Mott insulators, 
but it is recognized that the insulation performance, 
i.e. electron-to-electron correlative effect, increases 
along with the increase of a 2 . This is because the 
width of the effective Hubbard gap becomes larger as ct 2 
increases . 

Figs. 46, 47 and 48 show densities of states 
of half -filled electron systems (10001 electrons) in 
the nebula- like fractal-based complex structure under 
(ai, a 2 ) = (l, x) where x=0, 0.2, 0.4, 0.6 and 1. U=6 is 
used in Fig. 46, U=7 in Fig. 47, and U=8 in Fig. 48. 
In case of a 2 =l, the structure is the same as the 
simple fractal of a=l. In Fig. 46, for all values of 
a 2 , densities of states D(0) under co=0 are finite, and 
the electron system is in a metal phase. In Fig. 47, 
as a 2 increases, disappearance of D(0) occurs, that is, 
typical Mott transition occurs. In Fig. 48, for all 
values of ct 2 , the systems behave as Mott insulators. 
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but it is recognized that the insulation performance, 
i.e. electron-to-electron correlative effect, is 
enhanced as for all values of a 2 . This is because the 
width of the effective Hubbard gap becomes larger as a 2 
5 increases . 

Figs. 49, 50 and 51 show densities of states 
of half -filled electron systems (10001 electrons) in 
the nebula-like fractal-based complex structure under 
(ai, a 2 ) = (0.6, x) where x=0, 0.2, 0.4 and 0.6. U=6 is 
; s jl 10 used in Fig. 49, U=7 in Fig. 50, and U=8 in Fig. 51. 

m In case of a 2 =0.6, the structure is the same as the 

ro 

simple fractal of a=0.6. In Figs. 49 and 50, for all 
values of a 2 , densities of states D(0) under co=0 are 
*J finite, and the electron system is in a metal phase. 

15 As a 2 increases, D(0) gradually decreases, and the 

Q system approaches an insulator. In Fig. 51, for all 

ry 

□ values of a 2 , the systems behave as Mott insulators, 

rU 

but it is recognized that the insulation performance, 
i.e. electron-to-electron correlative effect, is 
20 enhanced as a 2 increases. This is because the width of 

the effective Hubbard gap becomes larger as a 2 
increases . 

Through the above analysis, it has been 
confirmed that various forms of correlative electron 
2 5 systems ranging from the metal phase to the insulator 

phase can be realized by changing (a lt a 2 ) . Especially 
when the condition of approximately D f <2.3 as the 
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fractal dimension in the initial stage of growth and 
the condition of approximately D f >2.7 as the fractal 
dimension of the final stage of growth are satisfied, 
remarkably excellent controllability can be obtained. 

Heretofore, the invention has been explained 
by way of some embodiments. The invention, however, is 
not limited to these embodiments but involves various 
changes or modifications within the technical concept 
and scope of the invention. 

As described above, according to the 
invention, by changing growth conditions with time 
during growth of a fractal structure to first grow it 
under growth conditions providing a low fractal 
dimension and thereafter grow it under different growth 
conditions providing a higher fractal dimension, it is 
possible to obtain a nebula-like fractal-based complex 
structure made up of a mixture of a plurality of 
regions different in fractal dimension from each other 
and to modulate and control the dimensionality of a 
material by a design method beyond those heretofore 
used to obtain conventional simple fractal properties. 
Moreover, the nature of various phase transitions 
occurring in these fractal structures can be 
controlled, for example, by adjusting the timing for 
changing the growth conditions. Furthermore, 
controllability can be improved by optimization of 
fractal dimensions. 



